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ABSTRACT 

We present main results of the investigation of the rare decay mode — > 77, in which the long 
distance contributions are expected to be dominant. Using the Heavy Quark Chiral Lagrangian 
we have considered the anomaly contribution which relates to the annihilation part of the weak 
Lagrangian and the one - loop vr, K diagrams. The loop contributions which are proportional to g 
and contain the ai Wilson coefficient are found to dominate the decay amplitude. The branching 
ratio is then calculated to be (1.0 it 0.5) x 10^^. Observation of an order of magnitude larger 
branching ratio could signal new physics. 
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In the past years the rare decays of B mesons came under the spothghts as the source of possible 
signals of new physics. In the meanwhile studies of rare D decays have received much less attention. 
Partially this is because theoretical investigations of D weak decays are rather difficult, also due 
to the presence of many resonances close to this energy region. The penguin effects on the other 
hand, which are very important in B and also in K decays, are usually suppressed in the case of 
charm mesons due to the presence of d, s, b quarks in the loop with the respective values of CKM 
elements. 

Nevertheless, D meson physics has produced some interesting results in the past year. Experi- 
mental results on time dependent decay rates of D'^ K^tt~ by CLEO and K^K^ and 

— > K~'K^ by FOCUS |^ have stimulated several studies on the — oscillations The 
recently measured D* decay width by CLEO [Q] has provided the long expected information on 
the value of D*D7r coupling. Among the rare D decays, the decays D — > V'j and D — > V{P)l^l^ 
are subjects of CLEO and FERMILAB searches On the theoretical side, these rare decays of 
charm mesons into light vector meson and photon or lepton pair have been considered lately by 
several authors (see, e.g., |11], for radiative leptonic D meson decay see [|l2|). The investigations 
of -D ^ Vj showed that certain branching ratios can be as large as 10~^, like for K*^^, 
Df /j^7 [^, However, the decays which are of some relevance to the — > 27 mode studied 
here, like p^^, — > ivy, are expected with branching ratios in the 10~^ range |13]. Thus, it 

is hard to believe that the branching ratio of the — > 27 decay mode can be as high as 10~^ in 
the Standard Model (SM), as found by 0]. Apart from this estimation, there has been no other 
detailed work on 27 prior to our analysis |T^. In addition to these theoretical studies there 

are experimental attempts to observe this decay rate done by CLEO and FOCUS 

On the other hand, in the B and K meson systems there are numerous studies of their decays 
to two photons. For example, the Bg — > 77 decay has been studied with various approaches within 
SM and beyond. In SM, the short distance (SD) contribution [|l^] leads to a branching ratio 
B{Bs 77) ^ 3.8 X 10~^. The QCD corrections enhance this rate to 5 x 10"'' Q. On the 
other hand, in some of the SM extensions the branching ratio can be considerably larger. The two 
Higgs doublet scenario, for example, could enhance this branching ratio by an order of magnitude 
|19|. Such "new physics" effects could at least in principle be dwarfed by long distance (LD) effects. 
However, existing calculations show that these are not larger than the SD contribution [20], which is 



typical of the situation in radiative B decays |21]. In the system the situation is rather different. 
Here, the SD contribution is too small to account for the observed rates of Ks — > 27, Kl — > 27 by 
factors of ~ 3 — 5 [22|, although it could be of relevance in the mechanism of CP-violation. Many 
detailed calculations of these processes have been performed over the years (see recent refs. [22| 
-[p5|] and refs. therein), especially using the chiral approach to account for the pole diagrams and 
the loops. These LD contributions lead to rates which are compatible with existing measurements. 

Motivated by the experimental efforts to observe rare D meson decays and noticing that 
Bs 77 offers possibility to observe physics beyond the SM, we undertook an investigation of the 
/)0 _j. decay [^]. Here we present only the main results of our analyses, while the details of 
our work are presented in |15|. 

The short distance contribution is expected to be rather small, as already encountered in the 
one photon decays [^, ^, hence the main contribution would come from long distance interactions. 
In order to treat the long distance contributions, we use the heavy quark effective theory combined 
with chiral perturbation theory (HQxPT) p^]. This approach was used before for treating D* 
strong and electromagnetic decays p7[1-p9[|. The leptonic and semileptonic decays of D meson were 
also treated within the same framework (see |27] and references therein). 

The approach of HQxPT introduces several coupling constants that have to be determined from 
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experiment. The recent measurement of the D* decay width Q has determined the 0*0-11 couphng, 
which is related to g, the basic strong couphng of the Lagrangian. There is more ambiguity, however, 
concerning the value of the anomalous electromagnetic coupling, which is responsible for the D*D^ 
decays (for further discussion on this point see ) . 

Let us address now some issues concerning the theoretical framework used in our treatment. 
For the weak vertex we used the factorization of weak currents with nonfactorizable contributions 
coming from chiral loops. The typical energy of intermediate pseudoscalar mesons is of order 
mi5/2, so that the chiral expansion p/A^ (for A-,^- ~ 1 GeV) is rather close to unity. Thus, for the 
decay under study we extend the possible range of applicability of the chiral expansion of HQxPT, 
compared to previous treatments like D* Dir, D* — > [^] or D* — > Djj ||2^, in which a 
heavy meson appears in the final state, making the use of chiral perturbation theory rather natural. 
The suitability of our undertaking here must be confronted with experiment, and possibly other 
theoretical approaches. 

At this point we also remark that the contribution of the order 0{p) does not exist in the 
1)0 _^ decay, and the amplitude starts with contribution of the order 0{p^). At this order 
the amplitude receives an annihilation type contribution proportional to the 02 Wilson coefficient, 
with the Wess-Zumino anomalous term coupling light pseudoscalars to two photons. However, 
the total amplitude is dominated by terms proportional to ai that contribute only through loops 
with Goldstone bosons. Loop contributions proportional to 02 vanish at this order. We point out 
that any other model which does not involve intermediate charged states cannot give this kind of 
contribution. Therefore, the chiral loops naturally include effects of intermediate meson exchange. 

The chiral loops of order 0{p^) are finite, as they are in the similar case of -fC — > 77 decays 



|22|-[^]. The next to leading terms might be almost of the same order of magnitude compared 



to the leading 0{p^) term, the expected suppression being approximately p^ / h^. The inclusion 
of next order terms in the chiral expansion is not straightforward in the present approach. We 
include, however, terms which contain the anomalous electromagnetic coupling, and appear as next 
to leading order terms in the chiral expansion, in view of their potentially large contribution (as 
in B*{D*) B{D)"fj decays considered in ||2^). As it turns out, these terms are suppressed 
compared to the leading loop effects, which at least partially justifies the use of HQ^PT for the 
decay under consideration. Contributions of the same order could arise from light resonances like 
p, K*, ao(980), /o(975). Such resonances are sometimes treated with hidden gauge symmetry (see, 
e.g., [^]), which is not compatible with chiral perturbation symmetry. Therefore, a consistent 
calculation of these terms is beyond our scheme and we disregard their possible effect. 

The invariant amplitude for 77 decay can be written using gauge and Lorentz invariance 

in the following form: 



M 



where M^") is a parity violating and M^^^ a parity conserving part of the amplitude, while ki(^2)^ 
ei(2) are respectively the four momenta and the polarization vectors of the outgoing photons. 

In the discussion of weak radiative decays q' — > 977 or q' q-y decays, usually the short (SD) 
and long distance (LD) contribution are separated. The SD contribution in these transitions is 
a result of the penguin-like transition, while the long distance contribution arises in particular 
pseudoscalar meson decay as a result of the nonleptonic four quark weak Lagrangian, when the 
photon is emitted from the quark legs. Here we follow this classification. In the case of 6 — > 577 
decay it was noticed that without QCD corrections the rate T{b — > s77)/r(6 — > 57) is about 
10"'^. One expects that a similar effect will show up in the case of c — > M77 decays. Namely, 
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Figure 1: One loop diagrams, not containing beta-like terms, that give nonvanishing contributions 
to the — > 77 decay amplitude. Each sum of the amplitudes on diagrams in one row Mj = Mi^j 
is gauge invariant and finite. Numerical values are listed in Table |^. 




Figure 2: Anomalous contributions to 77 decay. The intermediate pseudoscalar mesons 

propagating from the weak vertex are 7r^,r/, r/'. 



according to the result of |30| the largest contribution to c ^ n77 amplitude would arise from the 
photon emitted either from c or n quark legs in the case of the penguin-like transition c — > uj. 
Without QCD corrections the branching ratio for c ^ M7 is rather suppressed, being of the order 
10"^'' 1^ ^. The QCD corrections |^] enhance it up to order of 10~^. 

In our approach we include the c — > U7 short distance contribution by using the Lagrangian 

C = -^V^,V:,C'/^f^F,,m, {ua'^^l{l + ^,)c), (2) 

where rric is a charm quark mass. In our analysis we follow ^] and we take Cy^^ = (-0.7 + 
2i) X 10-2. 

The main LD contribution will arise from the effective four quark nonleptonic AC = 1 weak 
Lagrangian given by 

^ = "Tf E VugV:g[ai{qT''c)iuT^q) + a2{uT''c){qT^q)] , (3) 

q=d,s 

where T^ = 7^(1 — 75), aj are effective Wilson coefficients [^], and Vq^q^ are CKM matrix elements. 
At this point it is worth pointing out that long distance interactions will contribute only if the 
S'[/(3) flavor symmetry is broken, i.e. if rus / m^. Namely, due to VudV*^^ ~ —VusV*g and = 
the contribution arising from the weak Lagrangian disappears in the case of exact SU (3) flavor 
symmetry. 
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Figure 3: The diagrams which give nonzero amphtudes with one /3-hke couphng (denoted by •). 



Going from quark to meson level effective Lagrangian one uses heavy quark symmetry for c- 
quark and chiral symmetry of light quarks to construct HQ^PT Lagrangian This is then used 
to calculate the — > 77 deacy width to one loop order. Leaving out the details of our calculation 



(see [15|), we discuss the final results. 

The decay width for the — > 77 decay can be obtained using the amplitude decomposition 
in (I) 



(4) 



The main contribution to the decay width arises from the diagrams presented on Figs. ||, |^. 
The calculated amplitudes depend on the number of input parameters, as mentioned in [|l5|. The 
coupling g is extracted from existing experimental data on D* — > Dn. Recently CLEO Collab- 
oration has obtained the first measurement of D*~^ decay width T{D*~^) = 96 it 4 it 22 keV 
by studying the D*~^ — > D^7r~^. Using the value of decay width together with branching ratio 
Br{D*+ D^Ti+) = (67.7 ± 0.5)% one immediately finds at tree level that g = 0.59 ± 0.08. The 
chiral corrections to this coupling were found to contribute about 10% [27, |2^. In order to obtain 
the a coupling, we use present experimental data on Ds leptonic decays (/d ~ /^j^ = a/y/rrvo). 
In our calculation we take a = 0.31 GeV^/^ |l|]. For the Wilson coefficients ai we take 1.26 and 
02 = —0.47 [p3]. We give here the numerical results for the one-loop amplitudes in Table 1. 





Ml ^ [xlO-i° GeV] 


m|+^ [x 10-10 GeV-i] 


Anom. 





-0.53 


SD 


-0.27 -0.81i 


-0.16 -0.47i 


1 


3.55 +9.36i 





2 


1.67 





3 


-0.54 +2.84i 







4.41 +11.39i 


-0.69 -0.47i 



Table 1: Table of the nonvanishing finite amplitudes. The amplitudes coming from the anomalous and short distance 



{C^i^') Lagrangians are presented. The finite and gauge invariant sums of one-loop amplitudes are listed in the next 
r(±)\ 'pjjg nuijibers 1, 2, 3 denote the row of diagrams on the Fig.|l} In the last line the sum 



three lines (M?=^' = 
of all amplitudes is given. 



In the determination of D* D77 and B* S77 a sizable contribution from /3-like elec- 
tromagnetic terms |15] has been found [^]. Therefore, we have to investigate their effect in the 
77 decay amplitude. The nonzero parity violating parts of the one loop diagrams containing 
j3 coupling are given on Fig. |3|, while numerical results are presented in Table 2. 



Diag. 


m\ ^[x 10-10 GeV] 




/3.1 





-2.69 


/3.2 





2.69 


/3.3 





2.11 


/3.4 


0.88 


-0.007 


/3.5 





0.51 


/3.6 


-2.88 


-0.52 




-2.00 


2.09 



Table 2: Table of nonzero contributions of the amplitudes coming from the diagrams with P coupling. In the last 
line the sums of the contributions are presented. We use (3 = 2.3 GeV~^, rric = 1.4 GeV. 



Using short distance contributions, the finite one loop diagrams and the anomaly parts of the 
amplitudes and with numerical values of the amplitudes as listed in Table ||, one obtains 

Br{D^ 77) = 1.0 X 10-*^. (5) 

This result is slightly changed when one takes into account the terms dependent on (3. The branch- 
ing ratio obtained when we sum all contributions is 

Br{D^ 77) = 0.95 x 10"^. (6) 

By varying (3 within 1 GeV"^ < /? < 5 GeV"^ and keeping g = 0.59 it 0.08, the branching ratio 
is changed by at most 10%. On the other hand, one has to keep in mind that the loop contributions 
involving (3 are not finite and have to be regulated. We have used MS scheme, with the divergent 
parts being absorbed by counterterms. The size of these is not known, so they might influence the 
error in our prediction of the branching ratio. Note also that changing a would affect the predicted 
branching ratio. For instance, if the chiral corrections do decrease the value of a by 30% this would 
decrease the predicted branching ratio down to 0.5 x 10~®. 

We have presented here the results of the detailed calculation of the decay amplitude — > 77, 
which includes short distance and long distance contributions, by making use of the theoretical 
tool of Heavy Quark Chiral Perturbation Theory Lagrangian. Within this framework, the leading 
contributions are found to arise from the charged vr and K mesons running in the chiral loops, 
and are of order 0{p^). These terms are finite and contribute only to the parity violating part 
of the amplitude. The inclusion of terms of higher order in the chiral expansion is unfortunately 
plagued with the uncertainty caused by the lack of knowledge of the counterterms. As to the parity 
conserving part of the decay, it is given by terms coming from the short distance contribution, the 
anomaly and from loop terms containing the beta coupling, the latter giving most of the amplitude. 
The size of this part of the amplitude is approximately one order of magnitude smaller than the 
parity violating amplitude, thus contributing less than 20% to the decay rate. Therefore, our 
calculation predicts that the D — > 27 decay is mostly a parity violating transition. 

In addition to the uncertainties we have mentioned, there is the question of the suitability of 
the chiral expansion for the energy involved in this process; the size of the uncertainty related to 



5 



this is difficult to estimate. Altogether, our estimate is that the total uncertainty is not larger than 
50%. Accordingly, we conclude that the predicted branching ratio is 

Sr(L>° ^ 77) = (1.0 ± 0.5) X 10"^ (7) 

The reasonability of this result can be deduced also from a comparison with the calculated decay 
rates for the — > p{uj)j, which are found to be expected with a branching ratio of approximately 

10-6 |,0,|il. 

We look forward to experimental attempts of detecting this decay. Our result suggests that the 
observation of I? —> 27 at a rate which is an order of magnitude larger than (0), could be a signal 
for the type of "new physics" , which leads to sizable enhancement [^] of the short-distance c — > ti7 
transition. 
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